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1. Introduction 


The Long Duration Exposure Facility (LDEF) satellite was recovered in January 1990 after 
almost six years in space. LDEF was well-instrumented with ionizing radiation 
dosimeters, including thermoluminescent dosimeters (TLD’s), plastic nuclear track 
detectors (PNTD’s), and a variety of metal foil samples for measuring nuclear activation 
products. 1 In addition, the induced radioactivity produced in various spacecraft 
components provides information on the radiation exposure. 2 Analysis of these LDEF 
data by several groups is in progress under coordination of the LDEF Ionizing Radiation 
Special Investigation Group. 3 

The extensive LDEF radiation measurements (Fig. 1) provide the type of radiation 
environments and effects data needed to evaluate and help resolve uncertainties in present 
ionizing radiation models and calculational methods (Fig. 2). The LDEF data are 
particularly important to improving models for addressing radiation issues associated with 
Space Station Freedom since LDEF had the same altitude range (= 350 - 500) and orbit 
inclination (28.5°) as planned for the Space Station. In conjunction with the LDEF data 
analysis, a calculational program has been established at the NASA MSFC to provide 
calculational support to aid in LDEF data interpretation and to utilize LDEF data for 
assessing the accuracy of current models. A summary of the calculational approach is 
given in Fig. 3. The present report describes some initial results from this LDEF 
calculational study. 

The purpose of the calculations reported here is to provide some initial results to aid in the 
LDEF data interpretation - namely, to obtain a general indication of: (a) the importance of 
different space radiation sources (trapped protons, galactic protons, albedo protons, and 
albedo neutrons), (b) the importance of secondary particles, and (c) the spatial dependence 
of the radiation environments and effects expected within the spacecraft. These are only 
scoping estimates because several important approximations have been made in this initial 
work — e.g., a one-dimensional (aluminum slab) model of the LDEF spacecraft is used, 
and the angular variation of the incident radiation (particularly the trapped proton 
anisotropy) is not accurately simulated. Subsequent calculations are planned which will 
remove these approximations. 

The calculational method (described in Sec. 2) uses the High Energy Transport Code 
(HETC) 4 to estimate the importance of different sources and secondary particles in terms 
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Fig. 1. Ionizing radiation measurements aboard the LDEF satellite. 
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Fig. 2. Significance of LDEF satellite data for validation of ionizing radiation models and predictive methods. 
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Fig. 3. Approach for calculationaJ program in support of LDEF satellite ionizing radiation assessments. 



of fluence, absorbed dose in tissue and silicon, and induced radioactivity as a function of 
depth in aluminum. Radioisotope production from aluminum and stainless steel is 
computed for varying aluminum shielding thicknesses. Results from these transport 
calculations are given in Sec. 3. 

As indicated earlier in Fig. 2, the LDEF spacecraft had a fixed orientation, which allows 
measurement of the trapped proton anisotropy in the South Atlantic Anomaly region 
because the anisotropy is not “averaged out” by spacecraft motion as is usually the case. A 
model describing this anisotropy has been developed recently by Watts, et al. at MSFC. 5 
LDEF data should enable a definitive test of the model, and detailed transport calculations 
using a 3-D spacecraft model and anisotropic trapped proton spectra are planned to compare 
with the LDEF data. However, the induced radioactivity very near the spacecraft surface 
should be relatively insensitive to spacecraft geometry, and some approximate 1-D 
calculations have been made of the anisotropy of near-surface activation and comparisons 
made with preliminary LDEF data. These preliminary anisotropy comparisons are given in 
Sec. 4. 





2. Calculational Method 


Radiation transport calculations have been performed to obtain scoping estimates of the 
depth dependence of fluence, dose, and induced radioactivity produced in the LDEF 
spacecraft due to ionizing radiation exposure. 

The radiation sources considered are trapped (V an Allen belt) protons, galactic protons, and 
“albedo” neutrons and protons emanating from the earth’s atmosphere due to cosmic ray 
bombardment. The source spectra used as input for the transport calculations are shown in 
Fig. 4, and the procedure used in estimating these spectra is described in Appendix A. 
(The trapped electron spectrum is shown in Fig. 4 for comparison only. Since the trapped 
electrons are of low energy and produce effects very near the spacecraft surface, they are 
not considered in the transport calculations here.) 

Spectra for the different sources were assumed incident isotropically on one side of a slab 
of aluminum 100 g/cm 2 in thickness. This is, of course, an important approximation, not 
only because it neglects the 3-D shielding effects of the spacecraft but also because the 
actual angular distribution of the incident radiation is very different for the different 
sources, as illustrated in Fig. 5. 

The transport calculations were carried out using the SAIC version 6 of the HETC code. 4 
This code uses Monte Carlo methods to obtain a detailed simulation of the radiation 
transport (Fig. 6). At each nuclear collision during the transport process, a calculation of 
particle transport inside the nucleus is performed using a high-energy intranuclear-cascade- 
evaporation (ICE) model 7 to obtain the multiplicity, direction, and energy of all secondary 
particles (Fig. 7). For low-energy neutron (< 20 MeV) transport the high-energy ICE 
model is not applicable, but various experimental data libraries and transport codes are 
available in this low-energy region. For the calculations here, the low energy neutron 
source computed by HETC is coupled to the MORSE Monte Carlo 8 code for low-energy 
neutron transport (Fig. 8). 

The main output obtained from the transport calculations is depth-dependent fluence 
spectra. These spectra are folded with the response functions given in Appendix B to 
estimate the absorbed dose in tissue and silicon as a function of aluminum shielding 
thickness. While HETC provides radionuclide production directly as a natural outcome of 
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Fig. 4. LDEF exposure to ionizing radiation. Shown are cummulative, orbit-average differential and 
integral fluence spectra over the 5.8 yr. duration of the LDEF mission. 
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Fig. 5. Illustration of the nonuniform angular variation of LDEF exposure to ionizing radiation. Trapped 
proton exposure occurs in the South Atlantic Anomaly region where the flux is highly anisotropic at LDEF 
altitudes, with protons confined mainly in planes perpendicular to magnetic field lines and with in -plane 
asymmetry due to the East-West effect. Galactic protons are blocked out from below by the shielding effect of 
the earth. The angular distribution of albedo neutrons emanating from the earth's atmosphere is also 
geometrically constrained due to earth shielding. 
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Fig. 6. Overview of the High Energy Transport Code (HETC). 
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Fig. 7. Illustration of the two-step intranuclear-cascade -evaporation model used in the HETC transport 
code for computing secondary particles from high-energy ("spallation”) nuclear collisions. 


9 














Input Data Data Preparation Codes Radiation Transport Codes Analysis Codes Output 





8. 

CO 

5 

V 

H 

5^ 

3 

« 

ft 

2 

i 

X 


0 

o 

a 

e 

a 

i 

CO 


j 


10 


Fig. 8. SAIC version of the HETC code system, which includes high-energy proton and neutron transport, low-energy neutron transport, and gamma-ray transport. 



















the ICE model calculation, the statistical accuracy is generally poor when the product 
nucleus mass is far removed from the target nuclear mass. Since large target-product mass 
differences are of interest in the present problem (e.g., 7 Be from Fe in stainless steel), we 
have used an alternate procedure in which the HETC (and MORSE) computed fluence 
spectra are folded with available activation cross sections to estimate radionuclide 
production. Radioisotope production from aluminum and stainless steel were calculated 
using the cross sections given in Appendix C. 
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3. Results of Transport Calculations 


Results from the transport calculations are arranged to show the contributions of different 
radiation sources and the contribution of secondary particles in terms of various effects. A 
summary of the results are given in this section, with additional results given in 
Appendixes D-F. The table below is a guide to the location of various results in the report 
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The spatial dependence of the results are in terms of the areal density depth in aluminum 
from 0 to 100 g/cm 2 . To roughly relate these thicknesses to LDEF, the spacecraft diameter 
is 32 g/cm 2 , and the length is 68 g/cm 2 . (This is based on an average density obtained 
from the overall dimensions of 14 ft. diameter x 30 ft. long, a spacecraft structure weight 
of 8,000 lb., and a weight of 13,400 lb. for the experiments 1 .) 

Fluence 

Figure 9 compares the proton and neutron fluences (over all energies) for all sources. For 
the trapped proton environment, the fluence from secondary neutrons exceed the proton 
fluence for penetration depths ~ 10 g/cm 2 . The magnitude, and spatial dependence, of 
secondary neutrons from galactic protons is comparable to the secondary neutrons from 
trapped protons. 

Dose 

Figure 10 compares the importance of different sources in terms of the absorbed dose in 
tissue and in silicon. The trapped proton source dominates for penetration depths <50 
g/cm 2 . The albedo sources contribute at most a few percent. Additional results for 
secondary particle contributions to the absorbed dose are given in Appendix E. 
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Fig. 9. Depth dependence of proton (primary and secondary proton contributions) and neutron fluences over all energies produced 
by trapped proton, galactic proton, albedo proton, and albedo neutron environments during the LDEF mission time. The different 
environments are all assumed incident isotropically on one side (0 depth) of an aluminum slab 100 g/cm 2 in thickness. 
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Fig. 10. Importance of various radiation sources in terms of absorbed dose in tissue and silicon for 
the total dose over the duration of the LDEF mission (top graph) and as a per cent of the total dose at 
each depth from all sources (bottom graph). 
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Activation 


Figure 11 compares the contribution of different sources to 22 Na and to 7 Be production 
fiom aluminum. The galactic source contribution exceeds the trapped source contribution 
for depths > 50 g/cm 2 for 22 Na production and 2 25 g/cm 2 for 7 Be production. The 
relative importance of the galactic source, which has a harder spectrum, is expected to be 
higher for the higher threshold activation products, which is consistent with these Be 
vs. 22 Na results. Figure 12 shows that for the trapped proton source and the case of Na 
production the secondary neutron contribution becomes important at depths 2 30 g/cm . 

Results from calculations of radioisotope production in stainless steel are given in 
Appendix F. 


( nuclei / g ) Production ( nude! / g ) 



2 

Depth in Aluminum (g / cm ) 



Fig. 11. Importance of different sources in terms of 22jsfa and ^Be production from aluminum. 
The production is normalized for the total lifetime of the LDEF mission. 
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Contribution to 22 Na Production from Aluminum (per cent) 22 Na Production ( nuclei / g ) 




Fig. 12. Contribution of secondary particles in producing 22jsf a from aluminum by incident 
trapped protons; top graph for production over LDEF mission, bottom graph as percentage of 
total production from primary and secondary particles by trapped protons. 




4. Approximate Estimate of Activation Anisotropy 
and Comparison with LDEF Data 

We consider here the directional dependence of induced radioactivity near the LDEF 
spacecraft surface due to anisotropy of the trapped proton exposure. 

The results for 22 Na production from aluminum shown previously (Figs. 11 and 12) 
based on a 1-D geometry (with irradiation from one side) indicate that the dominate 
production mode near the surface is from primary trapped protons. While for the actual 
spacecraft geometry there may be some contribution from radiation entering the “opposite 
side”, we neglect this contribution for now and assume a 1-D geometry model of 
effectively infinite thickness. 

We have used anisotropic trapped proton spectra based on the Watts, et al. 5 model (Fig. 
13) and modified the 1-D MSFC straight-ahead proton transport code of Burrell 11 to 
compute activation products (using the activation cross sections of Appendix C) and to 
estimate the anisotropy of 22 Na production at small depths in aluminum. The results 
(Fig. 14) show a West/East anisotropy ratio for ^Na production that varies from about 
a factor of 2 near the surface to a factor of 3.5 at 10 g/cm 2 depth. 

These results are compared in Fig. 15 with recent preliminary measurements by Harmon, et 
al. 10 for the 22 Na activation of LDEF aluminum clamp plates. (These plates are 
relatively thin, 1.29 g/cm 2 , and located on the surface of the spacecraft.) Based on an 
approximate fit we have made to the data points (Fig. 15), the measured West/East 
anisotropy is 1.8 compared to a calculated ratio of 2.0. In comparing absolute 
magnitudes, the calculations are higher than the measurements by about 30% for directions 
in the vicinity of East, and higher by about 50% for directions in the vicinity of West. 
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Anisotropy of Trapped Proton Spectra 



Fig. 13. (Top) Predicted anisotropy of trapped proton flux > 100 MeV , from Ref. 9, using the 
anisotopy model of Watts, et al 5 . (Bottom) Trapped proton spectra in a differential solid angle 
looking West and looking East and the West/East ratio of the two spectra. Both graphs for orbit- 
average spectra (28.5° inclination, 450 km. altitude, solar minimum). 
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Fig. 14. Approximate calculation of activation of aluminum from LDEF exposure to trapped protons in the 
South Atlantic Anomaly using a 1-D spacecraft model. The anisotropy of the trapped proton spectra is taken into 
account using the model of Watts, et al 5 . The top graph shows the depth dependence of the activation as a function 
of angle in the horizonal plane (90° from zenith) perpendicular to the spacecraft axis . The bottom graph shows the 
predicted increase in directionality with depth in terms of the activation on the West vs. East side of the spacecraft. 
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Fig. 15. Comparison of approximate calculations for 22 Na activation of LDEF aluminum clamp plates (used on outer 
surface of spacecraft to secure experimental trays) with the preliminary measurements of Harmon, et al. 10 The 
calculated activation is at the mid-depth (0.64 g/cm 2 ) of the clamp and based cm a I-D model of the spacecraft. The 
measured anisotropy (activation on West vs. East side of spacecraft) based on the aproximate data fit shown is 1.8; 
the calculated anisotropy at this depth is 2.0. 
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Appendix A 

Sources of LDEF Ionizing Radiation Exposure 

Trapped Protons (Omnidirectional) 

The omnidirectional trapped proton spectrum calculated by Watts 1 was used as the trapped 
proton source for the HETC transport calculations. These spectra are based on the 
AP8MIN and AP8MAX trapped proton environment models 2 and the IGRF 1965.0 80- 
term magnetic field model projected to 1964, the epoch of the proton models. The 
cumulative flux over the duration of the LDEF mission was estimated by Watts by 
performing orbit average calculations (28.5° inclination, circular) at altitudes of 258.5, 
255.0, 249.9, 230.0, and 172 nautical miles, which took place on mission days 0, 550, 
1450, 1950, and 2105. A linear variation is assumed between time points. At altitudes of 
230 and 170 nautical miles, the solar maximum model (AP8MAX) was used, with the solar 
minimum model (AP8MIN) used for other altitudes. 

The resulting omnidirectional, altitude-average differential and integral cumulative trapped 
proton flux spectra over the duration of the LDEF mission are shown in Fig. A-l. For the 
one-dimensional transport calculation, one-half of this fluence was assumed to be incident 
isotropically on one side of the slab of material. While isotropy is a reasonable 
compromise for use in a one-dimensional approximation, the actual angular distribution is, 
as shown in the Sec. 4, highly anisotropic. 

Trapped Electrons 

The trapped electrons are of such low energy that they contribute significantly to the dose 
only at small penetration depths (S 0.5 g/cm 2 ) (Ref. 1) and do not contribute at all to 
radionuclide production. Thus, transport calculations for trapped electrons have not been 
made here, but the trapped electron spectra are given in Fig. A-2 for comparison with 
trapped protons. These trapped electron spectra were computed by Watts 1 using the 
AE8MIN and AE8MAX trapped electron environment models 3,4 and the same mission- 
averaged method as for the trapped protons given above. 


A-l 


Galactic Protons 


For the galactic proton spectrum we start with the analytic fit given by Adams, et al. 5 for 
the exomagnetospheric, time-dependent spectrum, which at solar minimum and solar 
maximum reduces to 


F(E) = 10 m (E/117500) 3 

where m = 6.52 exp {- 0.8 (logioE) 2 } - 4.0 

a = - 2.2 { 1 - exp [-b (logio E) 2-75 ] ) 

b = 0.1 17 at solar minimum and b = 0.079 at solar maximum. Here F has units of protons 
m 2 steradian' 1 Me V 1 and E is in MeV. This fit to the galactic spectrum (multiplied by 
4 jc steradians, converted to cm' 2 , and multiplied by the LDEF mission duration of 2114 
days) is shown as the “exomagnetosphere” spectrum in Fig. A-3. 


To take into account the effect of geomagnetic shielding at the LDEF orbit, we have used 
the geomagnetic field “transmittance fraction” given in Adams, et al. 6 for 30° inclination at 
400 km altitude, which is based on the cosmic ray trajectory tracing calculations of Shea 
and Smart for effective geomagnetic cutoffs over a world-wide, longitude-latitude grid at 
400 km altitude and the orbit averaging method of Heinrich and Spill 7 . This fraction of the 
exoatmospheric galactic protons transmitted through the geomagnetic field is shown in Fig. 
A-4, and the result of applying this transmission factor to the exomagnetospheric spectrum 
gives the curves labeled “LDEF orbit” in Fig. A-3. 

Another factor influencing LDEF’s exposure to galactic protons is the shielding effect of 
the earth’s “shadow”. The solid angle occulation is 

A £) = 2ft { 1 - [(R e + h) 2 - R 2 ] 1/2 /(Re + h)} 

where Re is the earth’s radius (6371 km) and h is the orbit altitude. For an average LDEF 
altitude of about 450 km, A Q/4ft = 0.32. Thus, 32% of the 4ft solid angle is blocked by 
the earth, and the incident proton directions are within ±110° about the zenith direction. In 
the transport calculations, the incident galactic proton flux was assumed incident 
isotropically over ± 90° about the target surface normal. 
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The galactic proton integral and differential energy spectra over the LDEF mission duration 
are given in Fig. A-5. 

Albedo Protons * 

Secondary protons produced in the earth’s atmosphere by cosmic rays can escape upward 
as “splash albedo” and become trapped in the earth’s magnetic field when the proton energy 
is below the geomagnetic cutoff. These protons are guided by the field to impact with the 
atmosphere in the hemisphere opposite to their formation, providing a “re-entrant” albedo. 

The splash albedo spectrum has been measured by several balloon flights at different 
latitudes. In particular, Wenzel, et al. 8 and Pennypacker, et al 9 measured the albedo 
spectrum in the 4 MeV to 1 GeV energy range at about 4 g/cm 2 residual atmosphere over 
Palestine, TX (42° N geomagnetic latitude, 4.5 GV geomagnetic cutoff). Measurements of 
, the proton albedo by the Cosmos-721 satellite (polar orbit, 210-240 km) have been 
reported by Kuznetsov, et al. 10 For a 4.5 GV cutoff they find a similar spectral shape as 
for the balloon flights but a factor of 4 higher intensity, which Kuznetsov, et al. attribute as 
possibly due to the different angular distribution of albedo protons at satellite vs. balloon 
altitudes. 

For the splash albedo calculations in the present work we have used a fit to these satellite 
and balloon measurements, with the magnitude of the balloon data increased by a factor of 
4 and the reported measurements per steradian multiplied by 2k to obtain an 
omnidirectional flux. These data and the fit used are shown in Fig. A-6, with the fit being 

<|> = 0.001 13 exp (- 0.0095E), 10 < E < 1 15 MeV 

= 0.79 E' 161 , 115 <E <2000 MeV 

where <j> has units cm’ 2 s " 1 MeV " 1 . This differential flux multiplied by the LDEF 
mission duration, together with the corresponding integral fluence, is shown in Fig. A-7. 


* We wish to thank J. Adams, Naval Research Laboratory, for providing background material on albedo 
proton measurements. 
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Neutron Albedo 


Some of the neutrons produced in the earth’s atmosphere by cosmic-ray 
bombardment escape the top of the atmosphere to constitute a neutron albedo. Several 
measurements and calculations of the neutron albedo near the top of the atmosphere have 
been made -- e.g., Fig. A-8. The results of Fig. A-8 are for the upward moving flux at 45 
km altitude, 42° N geomagnetic longitude, and solar minimum. We have fit the calculated 
spectrum as: 

<j)(E) = 0.047E ' 2-88 , 

= 0.40 exp (- 0.97E), 

= 0.15 E" 134 , 

= 0.0086 exp(- 0.045E), 

+ 0.0021 exp(- 0.0085E), 10 < E < 200 

= 1.95 E' 1,61 , 200 <E £3000 

where <)> has units cm ' 2 s' 1 MeV' 1 andEisinMeV. 

The analytic fit of Fig. A-7 is scaled as follows to obtain an estimate of LDEF exposure to 
albedo neutrons. The maximum geomagnetic latitude reached by the LDEF orbit is = 
40°, which is approximately the latitude corresponding to the spectrum of Fig. A-7. From 
measurements of the 1-10 MeV albedo flux dependence on magnetic latitude, the variation 
of the albedo flux over LDEF orbits (ratio of flux at = 40° to flux at 7^ = 0°) is about a 
factor of 3, and the ratio of the maximum albedo flux to the 28° inclination orbit-average 
flux is estimated to be a factor of » 2. Thus, while a detailed orbit integration has not been 
carried out to obtain the average LDEF exposure to albedo neutrons, the A. m = 42° 
spectrum of Fig. A-7 is multiplied by 0.5 as an estimate of the orbit-average exposure. To 
take into account altitude differences, l/i 2 scaling is assumed and the 45 km spectrum of 
Fig. A-7 is multiplied by 0.88 to obtain the spectrum at 450 km, which is approximately 
the average LDEF altitude. Finally, the flux of Fig. A-7 is multiplied by the LDEF on-orbit 
time (21 14 days) to obtain the albedo neutron fluence over the duration of the LDEF 
mission. The product of these scale factors (8 x 10 7 ) times the analytic fit curve of Fig. A- 
7 gives the estimate used for LDEF exposure to albedo neutrons (Fig. A-9). 


10' 5 < E<0.1 
0.1 £ Ecl.O 
1.0 £ E < 10 
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At altitudes of about 450 km the albedo neutron directions are restricted within a cone half- 
angle of 70° about the zenith because the earth shields neutrons of other directions. Thus, 
in the transport calculations only neutron directions within ± 70° about the slab normal were 
allowed. 

Summary 

Table A-l below summarizes the energy range and normalization for the different sources 
used as input for the transport calculations. Also indicated is the angular distribution range 
assumed in computing the source spectra per unit solid angle. 


Table A-l. Source parameters used for transport calculations. 


Source 

Minimum 

Incident 

Energy 

Maximum 

Incident 

Energy 

Omnidirectional Integral 
Fluence above E ^ 
(cm ‘ 2 over LDEF Mission) 

Range of 
Angular 
Distribution 
(steradians) 

Trapped Protons 

15 MeV 

600 MeV 

4.3 X10 9 

4jc 

Galactic Protons 

3.2 GeV 

100 GeV 

2.8 x 10 7 

2 n 

Albedo Protons 

15 MeV 

3.5 GeV 

2.3 x 10 7 

4 n 

Albedo Neutrons 

1 keV 

3.0 GeV 

7.4 x 10 7 

1.3k 
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Proton Energy (MeV) 


Pig, A-l. LDEF exposure to trapped protons, averaged over LDEF altitudes and cummulative over the LDEF 
mission duration, calculated by Watts 1 using the AP8MAX and AP8MIN environment models.2 Shown 
here are omnidirectional spectra; as discussed in the text, the trapped proton spectra at the LDEF altitude and 
orbit inclination are actually highly anisotropic. 



Electron Energy (MeV) 

Pig, A-2. LDEF exposure to trapped electrons, averaged over LDEF altitudes and cummulative over the LDEF 
mission duration, calculated by Watts 1 using the AE8MIN and AEMAX environment models. 3 - 4 
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Fig. A-3. Galactic proton spectra in interplanetary space (from Adams, et al.5) and at LDEF orbit 
after attenuation by geomagnetic field. 
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Fig. A -4. Transmission factor for galactic proton penetration of geomagnetic field; averaged over 30° 
inclination, circular orbit at 400 km. altitude (adapted from Adams, et al.^ ). 
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Fig. A-5. Cummulative galactic proton spectra over the duration of the LDEF mission, 
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Pig. A-6. Measurements of the "splash” proton albedo spectrum at balloon^ 1 ® and satellite^ altitudes. 
(The balloon data have been multiplied by a factor of four to get agreement with the magnitude of the 
satellite data.) Also shown is a fit to the data used here as input for the transport calculations. 
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pig. A-7. Cummuladve albedo proton spectra over the duration of the LDEF mission. 






Fig. A-8. HETC code calculations of the neutron albedo spectra from cosmic-ray bombardment of the 
earth’s atmosphere. ** Shown are flux (<(>) and current (J) spectra for the upward moving ( 2k ) and 
omnidirectional (4k) neutrons at 5 ,gfcm 2 residual atmosphere, 42° N geomagnetic latitude, and solar 
minimum conditions. Also shown are data from balloon flight measurements by Preszler, et al. 12 and 
by White, et al. 13 . 
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Pig, A-9, Cummulative neutron albedo spectra over the duration of the LDEF mission. The points (open 
circles) are from HETC code calculations of the differential spectrum at the top of the atmosphere (from Fig. A 
8), which have been renormalized for the LDEF orbit 
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Appendix B 

Dose Response Functions 

To estimate depth dependent doses the flux spectra at various depths from the transport 
calculations were folded with dose response functions. The response functions used (Figs. 
B-l and B-2) are "surface doses" for protons or neutrons incident normally on a slab of 
tissue or silicon. 

For protons incident on tissue, the response function was estimated using the stopping 
power approximation of Burrell 1 below 60 MeV and the transport calculation results of 
Zerby and Kinney 2 in the range from 60 to 400 MeV, Alsmiller, et al. 3 from 400 MeV to 
3 GeV, and Armstrong and Chandler 4 from 3 GeV to 100 GeV. For neutrons incident on 
tissue, results from Irving, et al 5 were used below 60 MeV, from Alsmiller, et al. 3 
between 60 MeV and 3 GeV, and from Armstrong and Chandler 4 between 3 GeV and 100 
GeV. 

For protons incident on silicon, the Burrell 1 stopping power approximation was used 
below 200 MeV with the HETC code kerma factor calculations of Zazula, et al. 6 used 
above 200 MeV. The response to neutrons is based on the DLC-31 data library 7 below 20 
MeV and the Zazula, et al. 6 calculations at higher energies. 
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Fig. B-l. Flux-to-dose conversion factors used for absorbed dose in tissue. 
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Appendix C 

Activation Cross Sections 


Michel and co-workers (e.g., Ref. 1) have developed a set of activation cross sections for 
neutrons and protons incident on various elements by using a combination of experimental 
data, semi-empirical methods, and nuclear models. Predictions of the spatial dependence 
of radioisotope production in thick composite targets using the Michel, et al. cross section 
set folded with flux spectra calculated by the HETC transport code are in very good 
agreement with experimental data for high-energy proton irradiations 2 . Thus, in these 
initial calculations we have used the Michel, et al. 2 activation cross sections shown in Figs. 
C-l through C-3. 
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Fig. C-2. Cross sections for ^Co, 54Mn and ^Be production from iron; solid line by protons, 
dashed line by neutrons. 











Appendix D 

Additional Fluence Results 


Figs. D-l through D-4 compare the depth dependent fluences (over all energies) for 
primary and secondary particles for LDEF exposure to trapped proton, galactic proton, 
albedo proton, and albedo neutron sources. . 

Fluence spectra of protons and neutrons from all sources are compared for 10 and 50 g/cm 
aluminum shielding depths in Figs. D-5 and D-6. 
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Fig. D- 1 . Secondary particle production for trapped proton source. 
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Fig. D-4. Secondary particle production for albedo neutron source. 







Fig. D-5. Comparison of fluence spectra at a depth of 10 g/cm 2 in aluminum from LDEF exposure to different ionizing radiation sources. 
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Fig. D-6. Comparison of fluence spectra at a depth of 50 g/ctrfl in aluminum from LDEF exposure to different ionizing radiation sources. 





Appendix E 


Additional Dose Results 

Figs. D-l through E-4 compare primary vs. secondary particle contributions to the 
absorbed dose in tissue for each of the LDEF exposure sources considered. 
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Fig. E-l. Secondary particle contribution to absorbed dose in tissue for trapped proton source 
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Fig. E-2. Secondary particle contribution to absorbed dose in tissue for galactic proton source. 
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Fig. E-3. Secondary particle contribution to absorbed dose in tissue for albedo proton source. 



Fig. E-4. Secondary particle contribution to absorbed dose in tissue for albedo neutron source. 
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Appendix F 

Results for Radioisotope Production from Stainless Steel 


Induced radioactivity measurements are being made for several LDEF components which 
are made of stainless steel, and calculated results are given here for several radioisotopes 
produced in thin stainless steel samples behind varying thicknesses of aluminum shielding. 
The stainless steel composition used in the calculations (75.3% Fe, 15.4% Cr, and 4.3% 
Ni, by weight) is based on post-flight x-ray fluorescence measurements (made at MSFC, 
and provided by A. Harmon, MSFC/SSL) of segments of the LDEF trunion. 

Figures F- 1 through F-3 compare the contributions from different sources to each of the 
radioisotopes considered. Figures F-4 through F-8 compare the primary vs. secondary 
particle contributions to the production of each isotope for trapped proton and galactic 
sources. 
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Fig. F-2. Contribution of various space radiation sources to the production of 5 ^Co (top graph) and 
54 Mn (bottom graph) from stainless steel, normalized for LDEF mission duration. 
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Fig. F-3. Contribution of various space radiation sources to the production of ?Be from stainless 
steel, normalized for LDEF mission duration. 
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Fig. F-4. Secondary particle contribution to 5^Co production from stainless steel by trapped 
proton (top) and galactic proton (bottom) sources, normalized for LDEF mission duration. 
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Fig. F-6, Secondary particle contribution to production from stainless steel by trapped 
proton (top) and galactic proton (bottom) sources, normalized for LDEF mission duration. 
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